November 1994 

TPR-94-34 

ADP-94-22/T162 



Spin-dependent nuclear structure functions: general approach 
with application to the Deuteron * 

S.A.Kulagin a \ W.Melnitchouk a > fc , G.Piller c and W.Weise a ' fe 
a Institut fur Theoretische Physik, Universitat Regensburg, D-930^0 Regensburg, Germany. 
b Physik Department, Technische Universitat Miinchen, D-85747 Garching, Germany. 
Department of Physics and Mathematical Physics, University of Adelaide, S.A. 5005, Australia. 



Abstract 



We study deep-inelastic scattering from polarized nuclei within a covariant 
framework. A clear connection is established between relativistic and non- 
relativistic limits, which enables a rigorous derivation of convolution formulae 
for the spin-dependent nuclear structure functions gf and g^ in terms of off- 
mass-shell extrapolations of polarized nucleon structure functions, g^ and 
g2 ■ Approximate expressions for g^ 2 are obtained by expanding the off- 
shell gi2 about their on-shell limits. As an application of the formalism we 
consider nuclear effects in the deuteron, knowledge of which is necessary to 
obtain accurate information on the spin-dependent structure functions of the 
neutron. 
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I. INTRODUCTION 



Polarized deep-inelastic scattering (DIS) experiments have in recent years yielded a num- 
ber of important and sometimes unexpected results. The measurements by the European 
Muon Collaboration (EMC) of the g\ structure function of the proton over a large range 
of values of the Bjorken scaling variable x, when combined with flavor non- singlet matrix 
elements from weak decays, provided information on the singlet axial charge of the proton. 
The small size of this resulted in the so-called proton "spin crisis", which prompted a seri- 
ous reanalysis of the very ideas behind the quark model and the simple parton picture of 
DIS. More recent experiments on the proton by the Spin Muon Collaboration (SMC) ]2[ 
and the SLAC E143 Collaboration || have allowed more refined analyses of the x and Q 2 
dependence of the g\ structure function P]-p|]. (For a recent review see Ref. J7J.) 

As an extra source of information, it is important also to measure the neutron polarized 
structure function, g™. Besides revealing the spin structure of the neutron itself, measure- 
ment of g™ is essential for testing the fundamental Bjorken sum rule. The absence of free 
neutron targets means, however, that light nuclei have to be used instead for this purpose. 
The SLAC E142 Collaboration M has in fact measured the structure function of 3 He, which, 
because of the preferential antiparallel polarization of protons in the 3 He nucleus, is believed 
to be approximately equal to the polarized structure function of the neutron. In addition, 
the SMC has recently measured the g\ structure function of the deuteron — combined 
with either the proton or neutron ( 3 He) data, this can be used as a valuable consistency 
check on the other measurements. 

To obtain accurate information on nucleon structure functions from nuclear DIS data, 
it is of course essential to reliably subtract any nuclear effects in the extraction procedure. 
Away from the shadowing region at small Bjorken x (x < 0.1), the standard method for 
investigating nuclear effects is the so-called convolution model, in which the nuclear struc- 
ture function is expressed as a (one-dimensional) convolution of the spin- dependent nucleon 
structure function and the nucleon momentum distribution in the nucleus. The convolu- 
tion model follows from the impulse approximation, Fig.[l], if one assumes that factorization 
between photon-nucleon and nucleon-nucleus scattering amplitudes translates into factor- 
ization between structure functions, although in a covariant framework this assumption is 
generally not justified. 



For He targets, nuclear effects have been investigated in Refs. [10,11 1, and for the 



deuteron in Refs. flQ , [12| -4T7f . In Refs. |lqjl7|| relativistic effects in the deuteron were also 
included, although still within the confines of the convolution model. It was shown in Ref. 
|T6| , however, that relativistic corrections necessarily lead to a breakdown of convolution for 
g^ when the full off-mass-shell structure of bound nucleons is incorporated. Although the 
convolution-breaking effects are not large (typically ~ 0.5% in the deuteron [0), in any 
self-consistent calculation they must be included. 

Aside from the relativistic complications, the situation is not completely clear even in 
the non-relativistic approaches. There exists in the literature [T^-|T4] a variety of results 
for convolution formulae for g±, the derivation of which is often based on early convolution 
models for unpolarized scattering ||18|1 , in which the issue of off-shell effects was not seriously 
addressed. The need exists, therefore, to derive convolution formulae for spin-dependent 
structure functions systematically in the non-relativistic limit. 
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In this paper we present an analysis of the polarized structure functions of nuclei, starting 
from a covariant framework, and working consistently to order p 2 /M 2 in the bound nucleon 
momentum. We demonstrate that in this limit one does indeed recover (two-dimensional) 
convolution formulae, although with different "flux factors" (polarized nucleon momentum 
distributions) compared to those found in the literature. Our formalism enables us to con- 
sider both the g± and g 2 structure functions (or equivalently the transverse structure function 
gT = gi + g 2 ) on a similar footing. For the latter, we find that g 2 receives contributions 
from the g± as well as from the g 2 structure functions of the nucleon. All of the formal 
results are valid in the Bjorken limit for spin 1/2 and spin 1 nuclei. 

Taking advantage of the weak binding of nucleons in nuclei, we also derive expansion 
formulae for gf T in terms of derivatives of g^ T . We concentrate on the specific case of 
the deuteron, where we present a detailed comparison of the expansion formula results with 
those of the full non-relativistic convolution, as well as with previous relativistic calculations. 
In addition, we estimate for the first time the nuclear effects that one needs to account for 
when separating the nucleon gj< (or g 2 ) structure function from deuteron data. This will 
be important in view of upcoming experiments which will for the first time measure g 2 for 
deuterium targets. 

The contents of this paper are laid out as follows: in Section |TJ we present the general 
covariant framework in which g^ T of a nucleus A are expressed in terms of the off-shell 
nucleon propagator and the virtual nucleon hadronic tensor; in Section [HI] we describe how 
the relativistic expressions can be reduced by making a non-relativistic expansion of the 
nucleon propagator in medium; Section [IV] deals with the details of approximate expansion 
formulae which are obtained by Taylor-expanding the off-shell nucleon structure function 
about its on-shell limit; application of the formalism to the case of deuterium is presented 
in Section |V|; finally concluding remarks are made in Section |VJ. 



II. RELATIVISTIC FRAMEWORK 
A. Definitions 

Inclusive deep-inelastic scattering of leptons from hadrons is described by the hadronic 
tensor 

W yv {P,q,S) = ±- f^e^(P,S\[J^),J v (0)]\P,S), (1) 

where P and q are the four-momenta of the target and photon, respectively, and the vector 
S is orthogonal to the target momentum, P-S = 0, and normalized such that S 2 = — 1. 
For spin-1/2 targets, S is simply the target polarization vector, while for the spin-1 case 
S is defined in terms of polarization vectors e™ such that S a (m) = —ie a (e m *,e m ,P) /Mt, 
where m = 0, ±1 is the spin projection along the axis of quantization, denotes the target 
mass, and we define e a (a, b, c) = e al3 ^ u a/3b^c u . The hadronic tensor can be decomposed into 
symmetric (s) and antisymmetric (a) parts, 

W^(P, q, S) = W$(P, q,S) + i W$(P, q, S). (2) 



3 



With unpolarized (charged) lepton beams one is sensitive only to the symmetric part, which, 
at leading twist, depends on the spin-independent F\ and F 2 structure functions, and, for 
spin-1 targets, also on the structure functions &i )2 when the target alone is polarized |I9| . 
If the lepton and hadron are both polarized, only the antisymmetric component of W^ u 
is relevant. For either spin-1/2 or spin-1 targets this is expressed in terms of the two 
independent, dimensionless structure functions gi^'- 



W$(P,q,S) 



1 



2P-q 



G p (P,q, S) = 2M T 



SP(g 1+ g 2 )-P' 



P-q 



92 



(3) 
(4) 



In the Bjorken limit (Q 2 , P-q — > oo), in which we work throughout, both structure functions 
gi and g 2 exhibit scaling, i.e. up to logarithmic QCD corrections they depend only on the 
ratio Q 2 /2P ■ q. 



B. Nucleon Tensor and Structure Functions 

In our covariant analysis it will be useful to work with the off-shell nucleon tensor W^ u , 
which is defined through the imaginary part of the forward photon scattering amplitude from 
an off-shell nucleon. In terms of W^, the hadronic tensor of an on-shell nucleon (p 2 = M 2 
where p is the nucleon momentum) is: 



W»(p,q,s)=Tr 



(^+M)^±^lw flu (p } q) 



(5) 



where s is the nucleon spin vector, and M is the nucleon mass. The antisymmetric part of 
W^v is given by an expression similar to that in Eq.(^), namely: 

HAp, q) = w q a G p (p, <?)• (6) 

zp ■ q 

In analyzing the tensor structure of it is convenient to expand in terms of 
a complete set of Dirac matrices, {/, 7", a a/3 , 7^75, 75}. The various coefficients in this 
expansion must be constructed from the vectors p and q, and from the tensors g a /3 and e^ a p. 
Terms proportional to q 13 do not contribute to W M ^, and are therefore not considered. The 
requirements of parity, time-reversal invariance and hermiticity restrict further the number 
of possible terms. In particular, it follows from parity invariance that transforms as an 
axial vector and the coefficient of the scalar (/) term must be zero. The requirements of 
time reversal invariance and hermiticity rule out the vector (7°) and the pseudoscalar (75) 
terms as well. Finally we find that can be written in terms of six independent structures: 

I nip) nii) \ 
Q\V, q) = (jpt + —ij P P 1 5 + G (7 V75 

+ 77 Pa + q a H 5 cr pa + — p^ l5 a aX p a q x , (7) 

\ M p-q J M p ■ q 
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where the coefficient functions G^ are constructed to be scalar, dimensionless and real 
functions of q 2 , p ■ q and p 2 . 

Substituting into Eq. @ we can express the on-shell nucleon [] structure functions g^ 2 
in terms of the coefficients (see also Ref. ||16|| ): 

g" = G iq) + G M + G {ap) - G {apq) , (8a) 
= -Gfa) + G {aq) + G {apq \ (8b) 

with the functions evaluated here at their on-shell values. Note that the term pro- 
portional to ^75 vanishes when tracing with the projection operator in Eq.(|5|), so that 
only five terms out of a possible six in Eq.(^) contribute to the physical nucleon structure 
functions. The structure ^75 could, in general, give non-vanishing contribution to structure 
functions of nuclei. However, as we shall see in Section |T|, only the above five functions will 
be relevant in the non-relativistic limit. 



C. Nuclear Structure Functions 



Discussions of nuclear effects in deep-inelastic scattering are usually framed within the 
context of the impulse approximation for the nucleons, see Fig.0. Other possible nuclear 
effects which go beyond the impulse approximation are final state interactions between 
the recoiling nucleus and the debris of the struck nucleon |[2Q|| , corrections due to mesonic 
exchange currents fl2lH23|, and nuclear shadowing 0. One may argue that complications 
due to meson exchange currents are less important here than in unpolarized scattering 
since their main contribution comes from pion exchange. Because it has spin zero, direct 
scattering from a pion constituent of a nucleus gives no contribution to spin-dependent 
structure functions. Also coherent multiple scattering effects, which are known to lead to 
nuclear shadowing, should not be important for large values of the nucleon Bjorken scaling 
variable x = Q 2 /2Mq . This is evident if one recognizes that the characteristic time scale 
1/Mx of the DIS process is smaller than the typical average distance between bound nucleons 
in the nucleus for x > 0.1. Based on these observations we consider the diagram in Fig.[TJ as 
a basic approximation. In this case the nuclear hadronic tensor can be written: 

W*(P, q,S) = J [dp] Tr [A(p; P, S) W^(p, q)] , (9) 

with [dp] = d A p/ (27r) 4 i. The function A(p; P, S) is the nucleon propagator inside the nucleus 
with momentum P and polarization S, 



Ap; P, S) 



d 4 Ce ip< (P,S T(N(£)N(0) 



PS 



(10) 



Throughout we define to be the average of the proton and neutron structure functions: 
< 2 = (#2 + #, a )/2. 

2 Some potential problems associated with the use of the impulse approximation for the g 2 structure 
function have also been discussed in Ref. [^4|] in the context of relativistic light-front dynamics. 
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Here iV(£) is the nucleoli field operator, and W^ip, q) is the hadronic tensor of the off-mass- 
shell nucleon, given by Eqs.(^) an d (0). 

The expression for the nuclear tensor in Eq.@ is covariant and can be evaluated in any 
frame. It will be convenient, however, to work in the target rest frame, in which the target 
momentum is P = (M A ; 0), and the momentum transfer to the nucleus, q = (q ; 0±, — |q|), 
defines the z-axis. For the g\ structure function it is natural to choose the spin quantization 
axis such that the target is polarized in a direction longitudinal to the momentum transfer, 
S = S 1 " = (0; 0j_, 1). Taking the W-fy component in Eq.@, which is proportional to gf, we 
find: 

xa9?(xa) = ^tJ Tr [A"(P) + Gz)] , (Ha) 

where xa = Q 2 /2P ■ q and x' = Q 2 /2p ■ q are the Bjorken variables for the nucleus and 
bound (off-mass-shell) nucleon, respectively. The nucleon propagator in the nucleus at rest 
with polarization S" is denoted by A"(p) = A(p; S"). 

The transverse spin-dependent structure function of the nucleus, g§ = gf + g%, is ob- 
tained by choosing the target polarization in a direction perpendicular to the momentum 
transfer, S = S 1 - = (0; S^, 0), and taking the and components in Eq.(^): 

XAsftixA) = ^fj [M x' Tr [AHP) G±] , (lib) 

where A ± (p) = A(p; S 1 - 1 ), and Q± represents the transverse spatial components of corre- 
sponding to the transverse quantization axis. From Eqs.( |lla] ) and ( llbj ) one can reconstruct 



the g^ structure function by taking the difference between g£ and gf. 

We should stress that the treatment culminating in the results of Eqs.flTT]) has been 
fully relativistic, and exact within the impulse approximation. In the literature one usually 



encounters formulations in terms of simple convolution formulae [pl)|- [l5| , p!7|| , in which the 
nuclear structure functions are expressed as one- dimensional convolutions of the nucleon 
momentum distribution in the nucleus, and the (on-shell) structure functions of the nucleon. 
To obtain the simple convolution result one requires two conditions to be satisfied: firstly, 
that the traces Tr [AG] in Eqs. flTTD factorize into completely separate nuclear and nucleon 
parts, and secondly that the nucleon component be independent of p 2 . For the twist- 2 
structure function gj 4 , it was shown in Ref. |I6| that in a relativistic treatment the off-shell 



degrees of freedom associated with bound nucleons in fact violate both conditions, leading 
to a breakdown of the simple convolution picture. A similar breakdown must also occur for 
the g£ structure function, since its twist- 2 contribution contains a component proportional 
to gf (see Section |V A| below). 

Clearly, in a relativistic theory one needs to go beyond the simple convolution formula- 
tion. On the other hand, the advantage of the convolution model is its ease of application. 
In the study of DIS from nuclei, especially light nuclei where typical binding energies are 
small, it may in fact be quite sufficient to treat the nucleus as a non-relativistic system. It 



was shown in Ref. |25| that for unpolarized nuclear structure functions, in which both of 
the above criteria are also violated relativistically [ ]2T)| , |2"T| ] , one can in fact recover factorized 
expressions in the non-relativistic limit. It is possible to define "off-shell nucleon structure 
functions", with the correct on-shell limits, which lead to (two-dimensional) convolution 
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formulae. In the next section we perform a similar non-relativistic reduction of the rela- 
tivistic expressions in Eqs.flTTD to establish whether a similar factorization is attainable in 
spin-dependent processes. 



III. NUCLEAR STRUCTURE FUNCTIONS IN THE NON-RELATIVISTIC 

LIMIT 

Our basic assumption in the remainder of the paper is that the nucleus is a non-relativistic 
system, made up of weakly bound nucleons interacting via the exchange of mesonic fields. 
This necessarily involves neglecting antinucleon degrees of freedom, and corresponds to 
bound nucleons in the nucleus being slow, |p| -C M . With these assumptions we will de- 
rive from Eqs.flTTD simplified (convolution) expressions for the spin-dependent g^ T structure 
functions of weakly bound nuclei. 



A. Non-Relativistic Reduction 



Following the procedure outlined for example in Ref. ||28|| , we can derive the relation 
between the relativistic nucleon field operator N and the non-relativistic operator ip. A 
detailed discussion of the non-relativistic reduction of N is given in Appendix The 
essential result is that, up to order p 3 /M 3 corrections, the operators N and ijj are connected 
via: 

mr^e^(^), , 12 , 



2M 



where p = — zV, and the renormalization operator Z — 1 — p 2 /8M 2 guarantees baryon 
number conservation. This result is valid for a wide range of meson-nucleon interactions, in 
particular for interactions with scalar and vector mesons as well as for pseudovector coupling 
to pions. Furthermore, it is explicitly interaction-independent. For the pseudoscalar 7riV 
couplings discussed in Refs. p2| , |29| , however, one finds explicit interaction dependence in 
Eq . (|12D , although the pseudoscalar interaction is generally considered less reliable than the 
pseudovector model, which we restrict ourselves to in this paper (see Appendix [S]). 

Consider now the consequences of applying Eq. (|T2|) to the traces in Eqs.flTTD. We start 
by writing the nucleon propagator (|Trj) as: 

-A a p{p) = -i Jdte tpot (T (N a (p,t)Wp(p,0))) , (13) 



2M, 



where N a (p, t) = J d 3 r e~ tp r N a (r, t) is the nucleon operator in a mixed (p, ^-representation, 
and the brackets denote an average over the nuclear state, (■■■) = (A\ ■■■ \ A) / (A\A) . Writ- 
ing the four- momentum of the bound nucleon as p = (M + £;p), we can introduce the 
non-relativistic nucleon propagator ^4. NR , which is usually defined as: 

A N J(e,p) =-iJ ^(T(^(p,t)^,(p,0))), (14) 
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where a, a' are the non-relativistic, two-dimensional nucleon spinor indices. The relativistic 
and non-relativistic nucleon propagators can then be related by substituting Eq. ([T2|) into 
Eq.©: 



>p- (15) 
Here the operator U translates two-component spinors into four-component spinors: 



Zip 



cr-p 

2W 



(16) 



Equation (|15D can be used to reduce the traces of the relativistic propagator A with the 
various Dirac structures in to expressions involving the non-relativistic propagator *4.nr : 



1 



where 



2M A 
1 

2M A 
1 



Tr [7o75*40)] = tr [<S A NK (e, p) 
Tr [7j75^(p)] = tr Sj A NR (e, p) 



2M A 
1 

2iW4 



Tr [i^a Qj A(p)\ = tr 
Tr [z75ffij-4(p)] = tr 



CT ' P \ 4 



1 ' 2M 2 
M 



A NR (e, P) 



<S 



5j 



cr-p 

M : 



P^ 



2M 2 



<7 P 

^ + 2M 2 " P * 



(17a) 
(17b) 
(17c) 
(17d) 

(18a) 
(18b) 



and i,j denote spatial indices. The trace "tr" is taken with respect to the spin variable 
in two-component space, tr O = O aa . All corrections to Eqs. (17) are of order |p| 3 /M 3 
or higher. One notices that the operators Sqj have a structure similar to the time and 
space components of the spin four-vector (0; tr) boosted to a frame in which the nucleon has 
momentum p. 

From Eqs. ( |I7|) and (|18~D one can determine the traces relevant for gf T in Eqs. ([11]), namely: 



1 



2M A 



-Tr 



p 2 -M 2 
2M 2 



tr 



Al R (e, p) (So + S : 



2M A 



-Tr 



A L {p)Q A 



+ G*»> (l + + (l - 4^ I I tr 



2M 2 



2M 2 



+ (-(*> + G H (l - + G<™> (l + l^S- , I tr 



2M 2 



2M 2 



*4nr0 ; p) S± 



(19a) 



(19b) 



where <S + S z and <S± are given by dl8|) , and p 2 « M 2 + 2M(e — p 2 /2M) is the squared 
nucleon four-momentum (the e 2 term is dropped since it introduces corrections of order 
p 4 /M 4 ). The operator T 2 is given by: 

where is the transverse component of the nucleon three-momentum vector: p = (p±,p z ), 
and defines the transverse spin quantization axis, relative to the photon direction (not to 
be confused with the spin operator S±). Note that in the non-relativistic limit the structure 
G^) does not contribute to (|19|). 

An important observation which can be made from Eq. (|T9|) is that the nuclear structure 
functions are expressed in terms of only two combinations constructed from the in 
Eq.©, 

g^x^) = G^ + GW + (G**> - G™) (l + P ^f] , (21a) 



2M 2 

,JW) = G W + C<"» (l + ^f) + G"* (l - P ^f) , (21b) 



so that also 

.2 \f2\ i ,;2 

92 { X ,P 



<W) = -G«> + SW 1 - !LZ*L) + G<-> 1 + !LZ* 



2M 2 ; V 2M 2 

= ^£W) - gf(x,p 2 ). ' (21c) 

These can be considered as definitions of the polarized nucleon structure functions in the off- 
mass-shell region (c.f. Eqs.@ above) in the vicinity of p 2 ~ M 2 . Note that in the p 2 = M 2 
limit they reduce directly to the free nucleon structure functions defined in Eqs. (13). 



B. Two-Dimensional Convolution 



The definitions of the off-shell structure functions in Eqs. (|2"T|) can now be utilized in 
deriving convolution formulae for gf T . After substituting Eqs. ([19]) into (|TT|), we make use of 
the analytical properties of the nucleon propagator for the integration over the momentum p. 
Namely, we close the contour of integration in the upper half of the complex e (or p ) plane 
and pick the poles of .Anr(£) which correspond to the time ordering 6(—t) (*J)^(0)ip(t)^. Fi- 
nally, one obtains simplified versions of Eqs. ([□]) which relate the polarized nuclear structure 
functions directly to those of nucleons: 

~ded 3 p 



x 



xg£(x) 



(2vr) 4 
ded 3 p 



tr 



P"{e,p) (S + S, 



x'g?(x' 



P 



tr 



P 



x g T {x . 



p z ) + tr 



P 



To 



x g 2 [x , p , 



(22a) 
(22b) 



where here the nuclear structure functions are expressed as functions of the standard Bjorken 
variable x = Q 2 /2Mq = xaMa/M. For a general polarization state, the nuclear spectral 
function, V(e, p), is defined by: 
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V aal (e, P ) = Ys^nAvWnAv) 27C6 ( £ - E o( A ) + E n(A ~ 1, -P)) 



(23) 



where the summation is performed over the complete set of states with A — 1 nucleons. The 
functions ip n ,<r{p) — ((A — l) n , —p\i/jo-(0)\A) give the probability amplitude to find in the 
nuclear ground state a nucleon with polarization o and the remaining A — 1 nucleons in a 
state with total momentum — p (n labels all other quantum numbers). The non-relativistic 
energies of the target ground state and the A — 1 residual nucleons are denoted E Q (A) and 
E n (A — 1), respectively. Summing over polarizations a and a 1 , the spectral function V(e, p) 
is normalized to the number of nucleons A: 



^?tr[P( £)P )] =A. 



(24) 



Equations (^) can be written in a more familiar form as two-dimensional convolutions 
of the off-shell nucleon structure functions in Eqs.([H|) and nucleon momentum distribution 
functions D(y,p 2 ): 



Jd P 2 J^D 1 (yy)g^(^,p^, 



y 



D T (y,p 2 )g»[-,p 2 ) + D T2 (y,p 2 )g»[-,p 



2\ „N I 



X 



y 



v 



(25a) 
(25b) 



where y = (po + p z )/M is the fraction of the light-cone momentum of the nucleus carried by 
the interacting nucleon. The nucleon distribution functions are given by: 



Di(y, P 2 ) 
D T (y,P 2 ) 
D T 2(y,p 2 ) 



ded 3 p 
ded 3 p 
ded 3 p 



V\e,p)(S + S z ) 5[y 



M + e + p z 



M 



8[p z -^), (26a) 



V ± (e,p)S±\5[y 
V ± (e,p)f 2 \s(y- 



M + e + p z 
M 

M + e + p z 
M 



5 [p 2 - fJ 2 
5 (p 2 - f 



(26b) 
(26c) 



where /1 2 = M 2 + 2M(e — p 2 /2M). We observe that is expressed entirely in terms of g^, 
while receives contributions from g^ as well as from g 2 . Consequently, the g 2 structure 
function will receive contributions from g± in addition to g 2 . 



IV. EXPANSION FORMULAE 

The behavior of the nucleon distribution functions in Eqs.(PB|) is governed by the nuclear 
spectral function V(e,p). The region of importance for V(e, p) is |p| < pp, \e\ ^ p 2 F /2M, 
where pp is a characteristic momentum which determines the momentum distribution of 
nucleons in the nucleus. For heavy nuclei this is the Fermi-momentum, pp « 300 MeV. For 
a light nucleus, such as the deuteron, the analogous parameter can be determined from the 
average kinetic energy T as yMT « 140 MeV. Therefore the nucleon distribution functions 
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( p6|) are strongly peaked about the light-cone momentum fraction y = 1 and the on-mass- 
shell point p 2 = M 2 . 

This property of the distribution functions allows us to obtain approximate expressions 
for the nuclear structure functions in Eqs.(]25|). The p 2 dependence of the off-shell structure 
functions can be first approximated by expanding g^T^/ViP 2 ) i n a Taylor series around 
p 2 = M 2 : 



N I X 2 



dp 2 



(27) 



p 2 =M 2 



where gi T [x/y) = g\T{ x /y,P 2 — M 2 ) is the structure function of the (physical) on-mass- 
shell nucleon (see Eqs.(||) and fl2T|)). Expanding g^ T (x/y)/y in Eqs. (f25|) around y = 1, 
integrating the result term by term, and keeping terms up to order e/M and p 2 /M 2 , we 
obtain simple expansion formulae similar to those used in the analysis of unpolarized nuclear 
structure functions [ 30, 21, 25 1. For the structure function g^ we find: 



irf(x) - Cf g»( X ) + C? (xg-(x))' + cf s (xg? (x))" + C? 



A 

where the derivatives are taken with respect to x, and the coefficients are: 

Sq + S z 



p 2 =M 2 



cf 
cf 



M 2 



Pz + e 
M 



S + S z — £ 



m v 2m; / ' 

with <S>o + 5 2 obtained from Eqs. (]T8|) . Here the averaging (• ■ ■) denotes: 

7%,p) . 



A 7 (2vr) 4 

Analogously, the structure function g£ is given by: 

^(s)«4%r(*) + Cf (xg?{x))' + C {2) x(xg$(x] 



(28) 

(29a) 
(29b) 

(29c) 

(29d) 

(30) 



r (3) dg«(x,p 2 



d\np 2 



+ 



c!$g?(x) + eg (xg?(x))'. 



p 2 =M 2 

(31) 



The coefficients C^) are identical to C[ l> except for the replacements (0)" — > (0) J ~ and 
S + S z — > Sj_, and the coefficients C^l are: 

Ct2 — {%■ 



*(0 



Ls T2 



Pz_ 
M 



(32a) 
(32b) 
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We stress that the matrix elements are calculated for longitudinally polarized targets, 
while for Cj) and Cj) 2 the target is polarized in a direction transverse to the photon direction. 
Once and g§ are calculated, the difference can be taken and the approximate result for g 2 
obtained. Corrections to g± 2 evaluated via Eqs.(PB"D and ( |5T| ) are of order p 3 /M 3 or higher. 

In the derivation of Eqs.(|28|) and (|3l|) one neglects the lower limit in the ^/-integration in 
Eqs.(P3|), namely the condition x/y < 1. From Eqs.(|2T)D one can easily see that this condition 
gives practically no restriction on the integration region in ( p9| , |32"D for 1 — x > pp/M. For 
heavy nuclei, the expansion formulae can be used safely up to x ~ 0.7. For light nuclei such 
as the deuteron, which we consider in the next section, Eqs.fl2~8|) and ([31]) are expected to 
be reliable for x < 0.8. 



V. THE DEUTERON 



In this section we consider the application of the results of Sections |IJ and [TV] to the case 
of a deuterium nucleus. As mentioned in Section |, deep-inelastic scattering from polarized 
deuterons is an important means of obtaining information on the polarized neutron structure 
functions g™ 2 - Extracting information on g™ 2 from deuterium data is only meaningful, 
however, if one has a reliable method of subtracting the relevant nuclear effects. 

The several previous attempts to account for nuclear effects in the deuteron (for the g® 
structure function) have not always yielded consistent results. Some early attempts |12| were 
made within a time-ordered framework in the infinite momentum frame, which unfortunately 
in practice was problematic due to the lack of knowledge about deuteron wavefunctions 
in this frame. Subsequent analyses ]13|] utilized convolution formulae obtained in direct 
analogy with the convolution model for unpolarized scattering. Namely, a one-dimensional 
convolution formula was used with the same non-relativistic "flux factor", (1 +p z /M), as 
appears in the unpolarized deuteron F 2 d structure function |]2llj25| , |3Tlj3^| . Other attempts 



TM were based on an operator product expansion at the nucleon level |22|, however these 



led to different operators to those in Eq. (|22a|) . It is important, therefore, to clarify which 



operators, or "flux factors", are relevant for g± when all terms up to order p 2 /M 2 are 
consistently kept. 

In addition to g®, there is considerable practical value in understanding the nuclear 
effects on g^ (and hence g 2 ), to which little attention has been paid thus far. In view 
of upcoming experiments which will for the first time measure the deuteron g 2 structure 
function [33], an estimate of the relevant nuclear corrections to g^ is urgently needed. 

As a guide to evaluating the most efficient method for the nuclear data analysis, a 
comparison of the results for gf 2 calculated using the convolution (^5|) and expansion (p8||3TD 



formulae will indicate the reliability of the latter approach to the deuteron, which to date has 
not been explicitly tested. Essentially the only ingredient needed to evaluate the coefficients 
in Eqs.(^8|) and (|3l|) , as well as the traces in the distribution functions of Eqs. (|26|) , is the 
deuteron wavefunction. Before discussing the details of the deuteron case, however, we must 
first fix the nucleon inputs that will be used in the subsequent numerical calculations. 
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A. Nucleon Structure Function Input 



For the structure functions and their derivatives we shall rely on experimental results 
where appropriate, and use model input where data are not yet available. For the proton 
and neutron g\ structure functions we use the recent parametrization from Ref. |34| of the 
SLAC |3§, EMC and SMC § proton, the SLAC-E142 neutron (Helium-3) §, and SMC 
deuteron || data. As an illustration of the quality of the fit, we plot in Fig.^(a) the xg^ n 
structure functions at Q 2 = 4 GeV 2 , compared with the data. 

For the g% structure function the study of nuclear effects is more problematic, since 
this receives contributions from both twist- 2 and 3 operators, the latter of which contain 
quark-gluon interactions and also explicitly depend on quark masses. Following the standard 
decomposition of g 2 into the different twist components, we write: 



g N { WW ){x) 



(33) 



where the twist- 2 part is given by the Wandura-Wilczek relation ||36|| : 



N(WW) 
92 



(x) 



x y 



(l-8(l-x/y)) g?(y), 



(34) 



and satisfies the Burkhardt-Cottingham sum rule [[37]: 



dx g 2 



N, 



x 



0. 



(35) 



The twist-3 piece (g 2 ) of g 2 is at present not very well determined at all. There is 



disagreement even about the magnitude and sign of its moments |38|j39 

m 



Based on a co- 



variant parton model approach, Jackson, Roberts and Ross [HDJ argued in favor of a very 
small twist-3 component. Within bag models, on the other hand, one finds []3~8|j41|| quite a 
sizable g 2 contribution compared with the Wandura-Wilczek term. An additional problem 
is how to relate the structure function calculated in the bag (or some other) model, which 
one assumes to be applicable at some low resolution scale Q 2 = 0(Aq CD ), to that appro- 
priate to DIS experiments (Q 2 ^ 5 GeV 2 ). One prescription j|2| is to simply assume the 
validity of the Altarelli-Parisi evolution equations down to very low Q 2 . Even within this 
pragmatic approach, while the evolution of the twist-2 component follows that of g±, which 
is understood, for the g 2 piece only the N = 2 and 4 moments can be handled exactly. The 
solution adopted in [|4[] was to use an approximate solution of the evolution equations that 
was derived in the large- Nc limit. 



Since the only available data [||3] on the g 2 structure function cannot yet unambiguously 
discriminate between the various models of g 2 , we will estimate the size of the nuclear 
effects for several models. To cover the potential range of results for g 2 , we take g 2 = 0, as 
suggested by Ref. j40|, and also the bag model predictions, evolved to Q 2 ~ 5 GeV 2 using 
the prescription adopted in Ref. ETJ. As an indication of the possible variation of the total 



g 2 structure function with different twist-3 components, we plot in Fig.|2](b) g 2 for these two 

which will be relevant in the actual evaluation 
2"%|,|3"1"D. One can see that while the effect of the 



cases. Also shown is the sum, 
of the nuclear structure function, 



9? 



9? + 9? 



Eqs. 



twist-3 component is certainly not negligible, it does not alter drastically the overall shape 
and sign of g 2 and g^ . 
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From the two-dimensional convolution equations ( p5|) it is clear that a consistent de- 
scription of nuclear effects to order p 2 /M 2 requires modeling in addition the p 2 dependence 
of the off-shell nucleon structure functions (^(^P 2 )- To this order of accuracy this can be 
achieved by determining the slope with respect to p 2 at the on-mass-shell point, Eq.(^). 
One could, for example, formulate g± 2 m terms of relativistic quark-nucleon vertex functions 



as described in Refs. ||lq|26| , p7|^4| , ^5| , and calculate the derivative directly. Alternatively, to 
obtain a quick estimate of the overall order of magnitude of the off-shell effect, we can ex- 



tend the model of Ref . [25] , which is based on a dispersion representation of the unpolarized 
nucleon structure function, to the polarized case. 

Within the latter approach, in the impulse approximation for the quarks the g± structure 
function can be written: 



^max( x iP ) 9 /j 9 9 

ds I dk p(k , SjP , x), 



(36) 



where k^^x^p 2 ) = —xs/ (1 — x) + xp 2 is the kinematical maximum of the quark momentum 
squared k 2 , s = (p — k) 2 is the center-of-mass energy squared of the "spectator" quark 
system, and p is the quark spectral function extended to the nucleon off-mass-shell region. 
Following Ref. we assume that p has no explicit x-dependence, and that the spectrum 
in s can be approximated by that calculated for a single effective mass s, p oc 5 (s — s). At 
moderate Q 2 ~ 5 — 10 GeV 2 one finds typically s ~ 2 GeV 2 []25[| . Taking a factorized p 2 and 
k 2 dependence in p then gives as a product of two functions: 



9i(x,p 2 ) = ipijp 2 ) T{k 2 max {x,p 2 



(37) 



normalized such that tp(M 2 ) = 1 and hence T ik 2 max {x , M 2 y) = gf(x). The explicit p 2 - 
dependence in the function tp(p 2 ) is dynamical in origin, while the p 2 -dependence in the 
function Tild^^ that enters through the upper limit of the k 2 integration is purely kine- 
matical. After some algebra, the derivative of gi(x,p 2 ) with respect of p 2 can then be 
written: 



dgi(x,p 2 



dp 2 



til?) 



d(p(p 2 



p 2 =M 2 



dp 2 



p 2 =M 2 



x (l 



X 



M 2 (I 



x 



(38) 



To determine the slope of the (p(p 2 ) at p 2 = M 2 we assume that the first moment of the 
non-singlet part of gf is not renormalized off-shell (at least in the vicinity of p 2 ~ M 2 ): 



x dx dg?(x,p 2 ) 
dp 2 



0. 



(39) 



p 2 =M 2 



The motivation for this condition is that since the axial U(l) anomaly [EB] is absent in the 
non-singlet sector, in the chiral limit the axial charge of the nucleon is a conserved quantity. 
For the singlet component of g± little is known about how the effects of the axial anomaly 
extrapolate into the off-shell region, although one would not expect dramatic consequences 
as long as p 2 ~ M 2 . To satisfy Eq.(]3"9"|) with s = 2 GeV 2 one needs a slope of approximately 
dip{p 2 ) I dp 2 \ p2=M 2 ~ —0.16 GeV -2 . For larger values of s, s ~ 3 GeV 2 , the slope becomes 
~ -0.1 GeV -2 . 
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The extension of the above off-shell model to the g^ structure function is more prob- 
lematic. For example, there is no known justification for a normalization condition such as 
in Eq.fl39"|) to be valid for . Furthermore, there is little knowledge about how the higher 
twist correlation, or final state interaction, effects involving the (nucleon-quark) "spectator" 
system would modify the effective mass s. For these reasons we postpone for the time being 
a more detailed discussion about the off-shell dependence in g% . 

With these inputs for the nucleon structure function we can now proceed to evaluate 
numerically the structure functions of the deuteron. 



B. Convolution Results 



For the deuteron, the traces in Eqs. (|26|) can be expressed in terms of the deuteron 
wavefunctions \I/ m=+1 (p) with spin projection m = +1 along the axis of quantization (see 
Appendix |FJ) []: 

tr [V(e,p) 0] D =*k(p)(0<P> + 0<*>)* +1 (p) 2n5(e-e D + p 2 /2M), (40) 

where eo is the deuteron binding energy. The expectation values in Eq. ( [40]) with the opera- 
tors O = Sq + S z , <Sj_ and T 2 determine the deuteron distribution functions. From Eq. (|B5|) 
(for the axis of quantization parallel to the z-direction) we find: 



tfj'ftp) (S + S z )^>(p 



,(ID, 



4V 



1 



2M 2 I 



2 UW , , 

u H t= ( 1 
V2 



2 I 3 ~2 



Pz 



Ml 2M 



u 



w 
V2 



1 2 



(41a) 



where u(p) and w(p) are the S- and .D-state wavefunctions in momentum space, and p z = 
\p\Pz, with p z = cos 6 (8 is the angle between p and the z-axis). 

For a deuteron polarized in a direction transverse relative to the photon direction, 
S = S 1 , we have: 



An 2 



+w z 



1 



2M 2 I 



^(1-3(^.8- 



u 



V2 



(P± • s- 



(p± • s 1 

2M 2 



u 



w 
V2 



1 2 



(41b) 



4tt 2 (p±-S ± ) j 3 / p x , 

' - 1 p z w w — v 2' 

2 V Mr V 



IP 



M 



(41c) 



3 Note that the argument in the deuteron wavefunction is the relative nucleon momentum, which 
in the deuteron rest frame coincides with the single particle momentum p. 
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Substituting these results into Eqs.(|26|), we can now evaluate the g® T structure functions 
of the deuteron. In Fig.|3] we plot the ratio g® / gf as a function of x, for Q 2 = 4 GeV 2 , using 
wavefunctions from several different models of deuteron structure pTF|-[4"9f| . The differences 
between the curves at small x are due mainly to the different deuteron D-state probabilities, 
namely 5.8% for the Paris, 4.3% for the Bonn (full model) and 4.7% for the Buck/Gross 
(pseudo-vector coupling) models. In calculating the curves in Fig|| we have made the 
approximation g±(x,p 2 ) ~ g^{x, M 2 ). Within the model of Ref. |TB| it was found that this 
approximation is good to within 0.5% for x < 0.8. We shall illustrate the effects of the 
dgf I dp 2 term evaluated using the off-shell model of Section |V A| in the next section. 

The transverse structure function ratio, g® / 9t i ls shown in FigfJ, at Q 2 = 4 GeV 2 , 
calculated for the Paris wavefunction [|47j (solid curve) 0. The result turns out to be quite 
similar to the ratio of the g\ structure functions in FigfJ (dashed curve), the main difference 
being at large x, where the gx ratio rises above unity earlier. The faster rise is even more 
pronounced for the ratio of the twist-2 components of gr (dotted curve). At intermediate x 
values {x < 0.5) the g® I 9t ratio appears mostly independent of the model for g^ , so that 
one may reasonably safely extract information on glp (and hence g 2 ) from the transverse 
deuteron data. 



C. Expansion Results 



In the range of x at which most of the data are taken, namely x ^ 0.6, the arguments 
given in Section |TV| would suggest that the expansion formulae in Eqs.(|28]) and ( |3lD should 
be excellent approximations to the full convolution results. Since these would simplify 
the analysis of the deuterium data, one may then take advantage of the simple expansion 
approximations in this region of x. 

To test the reliability of the expansion approach, we must calculate the coefficients , 
Cj) and C^l for a deuteron target. Using the matrix element in Eq. ( |41a| ) we obtain explicit 
expressions for the coefficients in Eqs.([29|) of the various terms in the structure function 
g°{x): 



cf 

cf 



" \ Pd ) m 1 + 1a ( To + ^ir To2 _ To T2 ) ' 



3M\ T °-— To2 ~W T2 , 



P D )\e D \ + 2T -T 2 



(42a) 

(42b) 

(42c) 
(42d) 



4 Note that plotting the ratio of g 2 structure functions is not very instructive since g 2 changes 
sign — to obtain the nuclear effects on g 2 alone one can subtract the g\ component determined 
above. 
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where Pd = J dpw 2 {p) is the P-state probability in the deuteron, with p = |p| (not to 
be confused with the momentum four- vector p used above). Here T , T 2 and T 02 repre- 
sent the average nucleon kinetic energies associated with each component of the deuteron 
wavef unction: 

r »=i *h^' (43a) 

d P f^w 2 (p), (43b) 

T 02 = J dp — u{p) w{p). (43c) 

To illustrate the role of the various terms in the expansion we plot in Fig.H the zeroth 
order contribution, proportional to C\ , together with the higher order terms, scaled by 
a factor 100. Evident at large values of x (x ~ 0.8) is the role of the second derivative 
term, proportional to C[ 2 \ which gives the characteristic rise (due to Fermi motion) of the 
structure function ratio Q\ I Q\ Fig.p. The trough in the ratio at x ~ 0.6 

arises from the single differential term proportional to which is large and negative in 
this region. Although it does not give rise to any uniquely distinct features in the structure 
function ratio, it is clear that the off-shell component (dotted curve in Fig.[|) is of the same 
order of magnitude as the other higher order corrections, and must be included in any 
precision analysis of nuclear effects in the deuteron. 

In Fig|5]we compare the performance of the expansion formula for the ratio g± / gf with 
several other methods of computation, neglecting for the moment the off-shell contributions. 
The simplest approach (and the one used by the SMC in their analysis ||) is to use a 
constant depolarization factor, (1 — 3/2 Pd), which roughly corresponds to the first term in 
the Taylor series in Eq . (p8|) . In fact, the order p 2 /M 2 correction to (1 — 3/2 Pd) is of the 
order of 10%, i.e. an overall correction to the structure function ratio of ~ 0.5%. At present 
this is still smaller than the uncertainty in Pd between the different models. Compared 
with the convolution results, one sees that the expansion formula works remarkably well for 
x < 0.7, where the two results are almost indistinguishable. As expected, for x ^ 0.7 the 
expansion curve overshoots the convolution result, which is understood from the fact that 
in the expansion formula one is neglecting the lower limit of integration, namely x, of D(y) 
over y, and replacing it by zero. 

To examine the effect of the off-shell C[ 3 ^ term in Eq.(|28|) on the g\ ratio, we plot in 



Fig.|7| the ratio with dg^ /dp 2 = (dashed) and with the slope determined through Eq.([3~8|) 
(solid) with the mass parameter s = 2 GeV 2 . The overall effect on the shape of the ratio 
is quite minimal, and the trend follows the shape of the off-shell component in Fig|| Note 
that within this model the off-shell curve at large x (x > 0.7) approaches the on-shell limit, 
although in this region the expansion approximation itself is no longer accurate. Within the 
relativistic model of Ref. |TB| (dotted curve), the magnitude of the (negative) off-shell effects 



was found to increase rapidly beyond x ~ 0.8, which is a further reason why the expansion 
curves tend to be larger at very large x. In the intermediate- a; region, on the other hand, 
there is little to distinguish all of the curves for 0.2 < x < 0.7. 

For the transversely polarized function g®(x), using Eqs. (|41b| , f41c|) we obtain: 

4 0) = Cf , (44a) 
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c a». fl _3 PD )M + J_( ro + ^ ro2 _^ 2l , (44b) 





£ *' = 3s( I4 + f T --^ r 'l- (44C) 

c£ 3) = CT, (44d) 

4° 2 } = ~ (t - ^T 02 + -T 2 ) , (44e) 

3M V 10 5 7 ' V ; 

\- (T 2 - v/2T 02 ) . (44f ) 



T2 



5M 



The various components of xg® in Eq.(^l|) are shown in Fig|| As in Fig.|], the higher order 
corrections are scaled by 100. (As discussed in Section |V A| , we do not consider the off-shell 
term here.) While small, the corrections proportional to xg^ and its derivative are still 
of the same order of magnitude as the -dependent terms. 

Finally, the relevant quantity needed for the extraction of the transverse nucleon structure 
function from g^ is plotted in Fig.|9] (solid curves). The ratio g® / 9t comes out to be very 
similar to the one obtained from the full convolution model (dashed curves). In Fig|| curves 
(i) have ~g~2 = 0> while curves (ii) include the twist-3 component. As a reference point, 
the result with a constant depolarization factor, (1 — 3/2 P D ), is also shown (dotted line). 
These results indicate that the expansion formula (|3T|) is quite a good approximation to the 
convolution model (|25b ) over nearly the entire x-domain of current experiments. 



VI. CONCLUSION 

We have presented a formulation of spin-dependent deep-inelastic scattering from spin 
1/2 and 1 nuclear targets. Starting from a covariant framework we have derived non- 
relativistic convolution formulae for the nuclear gf and g?{= g^ + 9^) structure functions 
in terms of polarized nucleon distribution functions and off-mass-shell extrapolations of the 
nucleon structure functions g^ T . It is known that relativistically the factorization of nuclear 
and nucleon parts of the total structure function, which is necessary for convolution, does not 
hold. Our results, however, are self-consistent to order p 2 /M 2 in the nucleon momentum, 
and represent the first systematic derivation of convolution formulae for polarized structure 
functions of weakly bound nuclei. To this order, while gf can be expressed in terms of g^, 
we find that g£ receives contributions from both the g% and g± structure functions of the 
nucleon. 

We have further utilized the fact that the nucleon momentum distributions in non- 
relativistic nuclei are strongly peaked around the light-cone momentum fraction y ~ 1 and 
the on-mass-shell point p 2 ~ M 2 . Expanding the virtual nucleon structure functions about 
these points we obtained simple expansion formulae for g± T valid for x ^ 0.7. 

The performance of the convolution and expansion approaches was examined for the case 
of the deuteron, which has direct practical implications for the extraction of the free neutron 
structure function from deuterium data. For the g® structure function good agreement was 
found between the expansion approximation and the non-relativistic convolution for all x 
below ~ 0.6. Differences between the non-relativistic convolution and previous relativistic 
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calculations become noticeable only for x ^ 0.7. At smaller x the source of the largest 
uncertainty is the non-relativistic deuteron D-state probability. 

We have also investigated for the first time the nuclear effects relevant for the extraction 
of the neutron structure function gj, (or g%) from measurements of the transverse structure 
function of the deuteron g®. Qualitatively these were found to be similar to the nuclear 
effects for g±, and for x ^ 0.5 largely independent of the details of the twist-3 component 

"f <£■ 

We can conclude, therefore, that for both and g^ the non-relativistic expansion 
formula provides a simple and, within current experimental accuracy, reliable means of 
analyzing nuclear effects in the deuteron. In future high-precision experiments ^ ne 
role of relativistic corrections as well as corrections to the impulse approximation itself may 
be more significant. These experiments should provide us with valuable guidance as to the 
relevance of relativistic effects in light nuclei, and where the impulse approximation may 
break down. 
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APPENDIX A: NON-RELATIVISTIC REDUCTION OF THE NUCLEON FIELD 

OPERATOR 

To derive the relation in Eq.(|T2|) between the relativistic and non-relativistic field oper- 
ators, we first write the relativistic nucleon field operator N in terms of upper and lower 
components, if and Xi respectively: 

= (ad 

In the extreme non-relativistic limit p — > 0, x — > 0, and at zeroth order in |p|/M, the 
non-relativistic nucleon field is simply given by the "large" upper component ip. We require, 
however, a relation which is valid to order p 2 /M 2 . 

The equation of motion for a nucleon field in the presence of a "potential" V can be 
written: 

(i#-M)N = V, (A2) 

where 

V = S+Y+ij 5 P + ^ (A3) 

describes the interaction of a single nucleon with mesonic fields produced by the surrounding 
nucleons. The first two terms in Eq . (|A3|) describe the coupling of a nucleon to scalar and 
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vector meson fields, S = g s a and = g^oj^ respectively, while the last two correspond to 
pseudoscalar (PS) and pseudovector (PV) 7rA-couplings, P = g^ir and = (g' 7T /2M) d^n 
(for simplicity we ignore isospin). 

Written for the upper and lower components, the Dirac equation (|A2|) reads: 

(ido-S- Vo-tr- A) if - (n- a + iP - A ) x = 0, (A4a) 
(n • a-iP-A ) v - (id + 2M - S + V + a • A)x = 0, (A4b) 

where 11= p— V, and p = — zV is the momentum operator. In the non-relativistic limit the 
driving term in Eq.( A4b ) is 2M. Expanding \ as a series in inverse powers of M gives: 



/ _ ldo + S - Vp - a ■ A \ (n-a-tP-A )ip 
X \ 2M J 2M ' 1 ' 

where the leading term is of order |p|/M, while the next term is of order |p| 3 /M 3 . For weakly 
bound nuclei, such as 2 H or 3 He, the interaction V is of the same order as the kinetic energy, 
p 2 /2M. Furthermore, the kinetic and potential energies are almost equal in magnitude while 
opposite in sign f], so that one can treat the different parts of the interaction (S, V and 
er • A) as all being of the order p 2 /M 2 . To this order, it is sufficient therefore to keep only 
the first term in the first parentheses in Eq. (^3) ■ 

From the equations of motion for the mesonic field we observe that the ratio of the 
spatial (V) to time (Vo) components of the vector current is N^fN j N^§N ~ |p| / M, so 
that | V/Vq\ ~ |p|/M. Furthermore, since the time component of the axial vector interaction 
is given by the time derivative of the pion field, A = (g' n /2M) don, it is reasonable to assume 
that Aq ~ (p 2 /M 2 )P. Therefore to order p 2 /M 2 one can neglect V and Aq in Eq . (|A5|) , so 
that the lower component x can be written: 

x = Jm {(T ' v ~ iP) ^- (A6) 



Substituting Eq. (|A6|) into Eq. ( |A4a| ) leads then to the Pauli-Schrodinger equation for the 



two-component nucleon spinor: 

ih<P = (j^ + V, (A7) 

^nr = S + V + a • A + ^.(VP) + ^ , (A8) 

where Vnr is the non-relativistic analog of V. Recalling that we assume all parts of the 
interaction being of the same order as the kinetic energy, p 2 /2M, we observe from Eq.( [A8| ) 
that the PS term should be of order P ~ cr-p, and so must be kept in Eq.( 



5 Note however that in relativistic models of nuclear matter, such as the Walecka model, some 
parts of the interaction V can be large but of opposite sign (e.g. S and Vo), leading to a small 
overall V. 
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Note also that to order p 2 /M 2 the nucleon density N^N receives contributions from the 
lower component, Xi m which case tp cannot be identified with the properly normalized non- 



relativistic nucleon field ip. Following Ref. [KJ, we introduce a renormalization constant 
Z such that ip = Zip, with Z determined by the particle number (charge) conservation 
condition: 

f dh^if; = J d 3 r N^N . (A9) 
Inserting x m Eq. (|A6|) into (|A9|) we find, to order p 2 /M 2 : 

Z = 1 " (p 2 + ff'(VP) + P 2 ) , (AlO) 

Therefore the renormalization constant depends explicitly on the PS pion-nucleon inter- 
action. In the non-relativistic limit the PV and PS couplings result in identical P-wave 
pion-nucleon interactions, as seen from Eq.([A8|). However, the PS coupling also generates 



a strong S— wave ttN interaction (the term P 2 /2M in Eq.(AJ)). Considered alone, the PS 
term leads to incorrect ttN scattering lengths, and its contribution is only cancelled by the 
introduction of a non-linear aim coupling. In our model therefore we consider only the PV 
coupling, which does not lead to the spurious S— wave interaction. In this case the non- 
relativistic expression for the four- component nucleon spinor is given by Eq . (|12"D , and the 
renormalization constant is interaction-independent. 



APPENDIX B: DEUTERON IDENTITIES 

For completeness we present here some definitions and useful identities for the deuteron 
which are used in Section IV A. 



The deuteron wavefunction in momentum space is defined as: 

*m( P ) = ^ Up) - W(p) Xl , m , (Bl) 



where p — |p| and m is the projection of the deuteron spin on the axis of quantization, Xi,m 
is the spin 1 wavefunction of the two-nucleon system, and S*i2(p) is the tensor operator with 
P = p/IpI- We define the S- and D-state wavefunctions in momentum space as: 



V2% 2 u(p) = / dr pr j (pr) u(r), (B2a) 



o 



V2tt 2 w(p) = / dr pr j 2 (pr) w(r), (B2b) 



o 



where u(r) and w(r) are the standard wavefunctions in configuration space, and normalized 
such that: 

dp (u 2 (p) + w 2 {p)) = 1. (B3) 
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If S is a unit vector along the direction of the spin quantization, then 

\ X { m (<r (p) + <x (n) ) X i,m = mS, 



(B4) 



where cr^^ are SU(2) Pauli spin matrices acting on the proton and neutron wave function 
respectively. Using some simple relations from the SU(2) algebra we find: 

*l(p) (t (p) + ^ (n) ) <Mp) = m ^ S (u 2 + ^= - w 2 ^ + (p • S) «, (w - v^u) . (B5) 
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FIGURES 



FIG. 1. Deep inelastic scattering from a polarized nucleus in the impulse approximation. The 
momenta of the target nucleus (P), virtual nucleon (p) and photon (q) are marked, and S denotes 
the nuclear spin vector. 



FIG. 2. Nucleon structure function input used in the calculation of nuclear structure func- 
tions, for Q 2 = 4 GeV 2 : (a) parametrization ]34| of the proton fi"|,p|,|35f and neutron || xg\ 
data; (b) isoscalar nucleon structure functions xg^ (= (xg^ + xgi)/2) (solid), xg% (dashed), and 
x 9t ( = x 9\ + x 92^j (dotted). In (b) curves (1) contain the twist-2 component of g^ only, while 
in curves (2) g^ has in addition a twist-3 contribution based on the bag model calculation of Ref. 



FIG. 3. Ratio of the deuteron to nucleon g\ structure functions calculated via the convolution 
formula in Eq. (|25a| ), with the Paris |^7| (solid), Bonn [Q (dotted) and Gross [^9| (dashed) deuteron 
wavefunctions. The latter, which also include small P-state components, are renormalized so that 
the S- and D-state wavefunctions alone are normalized to unity. 

FIG. 4. Ratio of the transverse gx deuteron to nucleon structure functions within the convolu- 
tion model, Eq. (|25b| ), for the Paris wavefunction |4^| , with (solid) and without (dotted) the twist-3 
component of g^. For comparison the g\ /g\ ratio (dashed) from Fig.^ is also shown. 



FIG. 5. Contributions to the deuteron xg® structure function from various terms in the ex- 
pansion formula, Eq.(28): zeroth order term (solid), first derivative (dashed), second derivative 
(dot-dashed), nucleon off-shell contribution (dotted). The latter three higher order terms are 
scaled by a factor 100. 



FIG. 6. Deuteron to nucleon structure function ratio, g^/g^, using the expansion formula 
(solid), compared with the convolution result (dashed) from Fig.||, and with a constant depolar- 
ization factor (1 — 3/2 Pd) (dotted), with Pjj ~ 5.8% from the Paris wavefunction [47|. 



FIG. 7. Deuteron to nucleon g\ structure function ratio from the non-relativistic expansion 
formula with (solid) and without (dashed) the off-shell component in Eq.(28), for the deuteron 
wavefunction of Ref. [49]. Shown also is the result of the relativistic calculation of Ref. [16] (dotted). 



FIG. 8. Contributions to the xg® structure function in the expansion formula, Eq.(31): (a) 
xg^ components — zeroth order term (solid), first (dashed) and second (dot-dashed) derivatives; 
(b) xg2 (triple-dot-dashed) and (xg^)' (dashed) components, on the background of the zeroth 
order term (solid) in (a). 
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FIG. 9. Transverse deuteron to nucleon structure function ratio using the expansion formula 
(solid), Eq.(|3l|), compared with the convolution model result (dashed). Curves (i) have = 0, 
while curves (ii) include the twist-3 component. The dotted curve indicates the constant depolar- 
ization ratio, (1 — 3/2 Pd), with Pjj ~ 5.8% from the Paris potential [47|. 
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